Efficient calculation of the antiferromagnetic phase diagram of the 3D Hubbard model 
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The Dynamical Cluster Approximation with Betts clusters is used to calculate the antiferromag- 
netic phase diagram of the 3D Hubbard model at half filling. Betts clusters are a set of periodic 
clusters which best reflect the properties of the lattice in the thermodynamic limit and provide an 
optimal finite-size scaling as a function of cluster size. Using a systematic finite-size scaling as a 
function of cluster space-time dimensions, we calculate the antiferromagnetic phase diagram. Our 
results are qualitatively consistent with the results of Staudt et al. [Eur. Phys. J. B 17 411 (2000)], 
but require the use of much smaller clusters: 48 compared to 1000. 
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The accurate and efficient solution of lattice Hamil- 
tonians such as the Hubbard model is a long standing 
challenge in the theoretical condensed matter commu- 
nity. These lattice models are routinely solved on a finite 
periodic lattice, for example with Monte Carlo, and the 
calculated properties extrapolated to the infinite limit. 
Due to the numerical expense in solving these models 
for large lattices, it is imperative to choose lattices that 
are efficient for the estimation and extrapolation of the 
physical properties of interest. 

In this paper we use the Dynamical Cluster Approxi- 
mation (DC A) Q (for a review see 0) to explore the 
antiferromagnetic instability in the 3D Hubbard model at 
half filling, with 
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where c^J (creates) annihilates an electron with spin a 
on site i, rii a is the corresponding number operator, t the 
hopping amplitude between nearest neighbors and 
U the on-site Coulomb repulsion. We solve this model 
on a series of finite clusters chosen according to the cri- 
teria proposed by Betts et al.0,Q- We obtain converged 
results extrapolating from clusters of up to only 48 sites, 
which are in agreement with the calculations of Staudt et 
al.0, who used conventional cubic lattices of up to 1000 
sites and obtained the Neel temperature via the specific 
heat. 

To solve the Hamiltonian (JTJ we utilized the DC A 4J. 
For a 3D system the DCA maps the original lattice model 
onto a periodic cluster of size N c — L\ embedded in a self- 
consistent host. Thus, correlations up to a range £ < L c 
are treated directly, while the longer length scale physics 
is described at the mean-field level. With increasing clus- 
ter size, the DCA systematically interpolates between the 
single-site dynamical mean-field result and the exact re- 
sult, while remaining in the thermodynamic limit. We 
solve the cluster problem using Quantum Monte Carlo 



(QMC) 8]. At half-filling there is no QMC sign prob- 
lem; the only systematic error in the Monte Carlo is the 
time step error, which can be extrapolated away. 

In order to calculate the phase diagram of the sys- 
tem in the thermodynamic limit, we employ the scaling 
ansatz 

£( r N CA ) = L c where T$ CA is the Neel tempera- 
ture obtained from a DCA calculation with a cluster of 
linear cluster size L c . This form is justified if we envi- 
sion the lattice as perfectly tiled by a periodic array of 
non-overlapping clusters. This system becomes ordered 
when the antiferromagnetic correlations of the cluster 
reach the linear cluster size. According to this ansatz 
C(TD CA ) oc |T£ ca - T N |-" a L c , so that 
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where Tn is the true antiferromagnetic transition tem- 
perature in the thermodynamic limit. The exponent is 
well-known for the 3D Heisenberg model, where one finds 
v k 0.710- 

Betts et al. 0, E| systematically studied the 2D and 3D 
Heisenberg models on finite size clusters and developed 
a grading scheme to determine which clusters should be 
used in finite size simulations. The main qualification is 
the "perfection" of the near-neighbor shells: a measure of 
the completeness of each neighbor shell compared to the 
infinite lattice. A perfect finite size cluster has all neigh- 
bor shells up to the k-th shell complete, the k-th shell 
is incomplete, and all shells k+1 and higher are empty. 
The absolute deviation from this criteria is defined as the 
imperfection. I.e. if the cluster neighbor configuration 
is as described, except that the k-1 shell is missing one 
entry, the cluster imperfection is one. The second qualifi- 
cation is the cubicity 6], C — max[c\, c\ Vs )max[c.2, c^ 1 ), 
where ci = Z l / 2 l/d and C2 = 2 1 / 2 /// are defined by the 
geometric mean of the lengths of the four body diago- 
nals of the cluster, d = (didadsd^) 1 ' 4 , the six face diago- 
nals / = {hhhhhh) 1 ^, and the edges I = (hhh) 1/3 . 
C = 1 for a cube, so the deviation of C from one is a mea- 
sure of the cubic imperfection. In finite size scaling calcu- 



2 



lations of the order parameter and ground state energy, 
they found that the results for the most perfect clusters 
fall on a scaling curve, while the imperfect clusters gen- 
erally produce results off the curve. We generated addi- 
tional 3D clusters following these guidelines for clusters 
larger than the 27 site clusters previously published)^, 
from which we adopt the labeling conventions and cluster 
geometries. In Table[I]we list clusters up to 70 sites, their 
perfection and cubicity. In each case, we chose either the 
bipartite (labeled B) or non-bipartite cluster (labeled A) 
with the smallest imperfection and cubicity closest to 
one, in this order of priority. For example, the 38 site 
cluster 38B is bipartite, perfect, consisting of only com- 
plete neighbor shells, and has a cubicity of 1.087. Since 
we are interested in a calculation of T^ CA , we utilized 
only bipartite clusters in the present calculations. 

To obtain the antifcrromagnctic phase diagram we per- 
formed a series of DCA calculations as a function of U/t, 
cluster size, and Monte-Carlo time step At. For a given 
U/t and cluster size, we calculated T® CA by finding the 
divergence of the staggered susceptibility as a function 
of At, and extrapolated the value obtained to At = 0. 
As an example, we show in Fig^the Neel temperature 
T£ ca (At) for an 18 site cluster for U/t = 8. One finds a 
significant At dependence which makes an extrapolation 
to At = mandatory. 
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FIG. 1: Tjv versus At 2 when U/t = 8 for cluster 18A 

Performing this extrapolation for the series of bipartite 
clusters from Table [I] for U/t = 8, we obtain the values 
for Tj° CA collected in Fig. (full circles). For comparison 
we also included the results for a finite t ■ At = 1/4 
(open circles). This unextrapolated data actually lies 
above the Heisenberg result of T/t=0.48. One clearly 
sees that a proper scaling to At = is necessary to obtain 
both the correct qualitative and quantitative behavior of 
T$ GA (N C ). The full curves in Fig. were obtained with 
the scaling ansatz © using the v for the 3D Heisenberg 
model. It yields a linear scaling curve within our error 
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TABLE I: 3D cluster geometries, imperfection and cubicity of 
the best non-bipartite (A clusters) and bipartite (B clusters). 
The a,i denote the cluster lattice vectors. 
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FIG. 2: Cluster size scaling of Tn when U/t — 8 and i At = 
1/4 (open circles) and the result extrapolated to tAr = (full 
circles) as in Fig. 



To assess the value of Betts clusters, we also study 
two bad clusters, 16Z and 26Z, identified in Table [H] 
Although these clusters are bipartite, they are highly im- 
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Imperfection Cubicity 



16Z ( 2, 0, 0) ( 0, 2, 0) ( 0, 0, 4) 7 1.209 

26Z ( 1, 2, 3) ( 3, 3,-2) ( 3,-2, 3) 14 1.295 

TABLE II: 3D cluster geometries, imperfection and cubicity 
of two poor quality bipartite clusters. 

perfect. Both are missing independent neighbors in the 
first shell (each have 4; whereas a complete first shell 
has 6 neighbors). As a result of the periodic boundary 
conditions on the cluster, this causes the near-neighbor 
fluctuations to be over estimated. As a result, the esti- 
mates of Tn from these clusters, shown Fig.[2]for a finite 
t ■ At — 1/4 (open triangles) and for the data extrap- 
olated to At = (filled triangles), fall well below the 
scaling curve established by the best cluster geometries 
listed in table |U In general, in this and in other cal- 
culations, we find that the less perfect clusters tend to 
overestimate the effects of fluctuations. 

Finally, Fig [3] displays the calculated antiferromagnetic 
phase diagram obtained from the DCA and extrapolated 
to At = and N c = oo (open circles with error bars). 
For comparison, we included results from other methods: 
The dynamical mean-field approximation (DMFA, full 
circles), Staudt et al.Q (full curve), second order pertur- 
bation theory (SOPT, dotted curve) [Tol ITT| . the Heisen- 
berg model (dashed curve) [J^ an d the Weiss mean-field 
theory for the Heisenberg model (dash-dotted curve). We 
took J = At 2 /U for both Heisenberg calculations. The 
results from Staudt et al. are reproduced with good accu- 



racy, but with much smaller clusters. The DMFA result 
is obtained through the methods described above when 
N c = 1. Both the DMFA and the Weiss mean field are 
local approximations which neglect the effect of non-local 
fluctuations. As expected, they agree in the strong cou- 
pling regime, U > 12t = W (W is the bandwidth). Both 
DMFA and SOPT are only accurate at small U/t, indi- 
cating that non-local fluctuations are not important for 
small U. At large U/t the DCA results for Tn approach 
the curve for the Heisenberg model, as expected. How- 
ever, for intermediate and large values of U/W, the de- 
viation between the present results and the mean-field 
results is as large as 30% or more, indicating that the 
effects of non-local fluctuations are significant. 
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FIG. 3: Antiferromagnetic phase diagram of the 3D Hubbard 
model from our results and different approximations. 



In conclusion, we have calculated the antiferromagnetic 
phase diagram of the 3D Hubbard model at half filling us- 
ing the Dynamic Cluster Approximation and Betts clus- 
ters. Well converged results are found for relatively small 
cluster sizes due to the well-chosen geometries of these 
clusters. The dramatically increased efficiency of these 
clusters compared to typically used cluster geometries, 
such as cubic lattices, suggests that these clusters should 
be more widely used for lattice calculations. 
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